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ABSTRACT 

We report on two epochs of very-long-baseline interferometry (VLBI) observations of the Type Ib/c 
supernova, SN 2008D, which was associated with the X-ray outburst XRF 080109. At our first epoch, 
at t = 30 days after the explosion, we observed at 22 and 8.4 GHz, and at our second, at t = 133 days, 
at 8.4 and 5.0 GHz. The VLBI observations allow us to accurately measure the source's size and 
position at each epoch, and thus constrain its expansion velocity and proper motion. We find the 
source at best marginally resolved at both epochs, allowing us to place a 3c upper limit of ~ 0.75 c on 
the expansion velocity of a circular source. For an elongated source, our measurements are compatible 
with mildly relativistic expansion. However, our 3tr upper limit on the proper motion is 4 /ias d , 
corresponding to an apparent velocity of < 0.6 c, and is consistent with a stationary flux centroid. 
This limit rules out a relativistic jet such as a gamma-ray burst jet away from the line of sight, 
which would be expected to show apparent proper motion of > c. Taken together, our measurements 
argue against the presence of any long-lived relativistic outflow in SN 2008D. On the other hand, 
our measurements are consistent with the nonrelativistic expansion velocities of < 30, 000 km s^ 1 and 
small proper motions (< 500 km s _1 ) seen in typical supernovae. 

Subject headings: supernovae: individual (SN2008D) — radio continuum: general — gamma rays: 
bursts 



1. INTRODUCTION 

The X-ray outburst, XRF 080109, was serendipitously 
discovered during Swift X-ray obs ervations of the galaxy 
NGC 2770 (jSoderberg et al.ll2008f l. which is at a distance 
of 28 Mpc ( H a = 70 km s" 1 M pc" 1 ; z = 0.006494 from 
HyperLeda. iPaturel et al.ll2003l ). Shortly thereafter, op- 
tical spectroscopy revealed a Type Ib/c supernova, SN 
2008D, associated with the outburst (IBlondin et alJ2008t 
ISoderberg et alj 120081 : IValenti et alj l2008bl laD. Radio 
emissio n was detected using the NRAO 5 Very Large 
Array (|Soderberd 12008). and we undertook first epoch 
very-long-baselinc interferometry (VLBI) observations in 
2008 February, with preliminary results published in 
ISoderberg et al.l ((2008). We report here on our second 
epoch VLBI observations in 2008 May and a full analy- 
sis of both sets of VLBI observations. 

XRF 080109 was the first detection of such an X- 
ray outburst associated with a supernova (SN). It had 
a timescale of ~600 s, and has been interpreted both 
as the shock breakout emission from the SN (e.g., 
ISoderberg et al.ll2008t IChevalier fc FranssorJl2008ft. an d 
as an engine- driven event (e.g. JXu et al. 2008; Li 2008). 



It had a peak luminosity of ~ 6. 1 x 10 J erg s and a 
total energy of - 1.3 x 10 46 erg (jSoderberg et al.ll2008h . 
The X-ray spectrum seems to be somewhat better fitted 
with a power-law (photo n index = 2.3) than with a black- 
body (kT = 0.73 keV. ISoderberg et all 12001 : IXu et all 
l2lm lLl l2008l : lModiaz et al.ll2008h . 
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Although SN 2008D was first identified as a Type Ic, 
it su b sequently transitione d to Type lb (jModjaz et alJ 
120081: IValenti et all l2008afh Ear ly spectra exhibited 
fairly broad featur es (see, e.g., IBlondin et al.1 120081 : 
IValenti eT al. 2008b), although not as broad as seen in so- 
called "hypernovae" associated with gamma-ray bursts 
(GRBs) such as SN 1998bw and SN 2006aj. It has 
been suggested that there may have been r elat ivistic 
ejecta present in S N 2008D (iDado et al]l2008t Pl2008t 
iMazzali et ail 120081: IXu et al.ll2008ft . The high angular 
resolution of VLBI allows a direct measurement of the 
angular size, as well as the proper motion, and thus al- 
lows us to constrain the velocities present in the source. 

2. OBSERVATIONS 

On 2008 February 8, {t = 30 d, where t = is the 
explosion time, which was 2008 January 9), we observed 
SN 2008D using the NRAO Very Long Baseline Array 
(VLBA; 10 antennas, each of 25 m diameter), while on 
2008 May 21, we used the phased Very Large Array, 
(VLA; 130 m equivalent diameter), and the Effelsberg 
Radio Telescope 6 (100 m diameter) in addition to the 
VLBA. Further details of the observing runs are given in 
Table HJ The data were correlated with NRAO's VLBA 
processor, and the analysis carried out with NRAO's As- 
tronomical Image Processing System (AIPS). The ini- 
tial flux density calibration was done through measure- 
ments of the system temperature at each telescope, and 
then improved through self-calibration of the reference 
sources. 

For both epochs, we phase-referenced to 
VCS1 J0919+3324 (also JVAS J0919+3324 and 
IVS B0916+336; hereafter just J0919+3324, 
iBeaslev et al.1 |200"2I) . We used a cycle time of ~5 min, 

6 The 100 m telescope at Effelsberg is operated by the Max- 
Planck-Institut fur Radioastronomie in Bonn, Germany. 
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with ~3.7 min spent on SN 2008D, except for the obser- 
vations at 22 GHz, where we used a somewhat shorter 
cycle time of ~3 min with ~1.7 min on SN 2008D. 

3. RESULTS 

On 2008 Feb. 8, we detected SN 2008D/XRF 080109 at 

8.4 GHz, with a total flux density of ~1.5 mJy and an im- 
age rms background of 90 //Jy bm _1 . At 22 GHz it was 
not detected. The expected 22-GHz peak brightness is 
< 0.8 mJy bm _1 , depending on the degree of resolution. 
Our image had an rms brightness of 0.15 mJy bm _1 . 
However, since both positive and negative extrema were 
~ 0.65 mJy bm _1 , no peak could confidently be associ- 
ated with SN 2008D. 

On 2008 May 21, we detected SN 2008D/XRF 080109 
at both 8.4 and 5.0 GHz. Figure [1] shows an image of 
SN 2008D at 5.0 GHz, where the total flux density was 
320 /iJy and the rms background level was 21 /iJy bm _1 . 
At 8.4 GHz, the total flux density was 160 //Jy and the 
rms background level 22 //Jy bm . The considerably 
higher sensitivity of the May observations was due to the 
wider bandwidth, the inclusion of two large apertures, as 
well as the somewhat longer observing time. The source 
is very marginally resolved at our FWHM resolutions of 

2.5 x 0.8 mas at p. a. —15° at 5 GHz and 1.8 x 0.6 mas at 
p. a. —17° at 8.4 GHz. In particular, no jet-like extension 
is seen. 

3.1. Modelfitting to Determine Source Sizes 

In order to accurately estimate the angular size, we 
turned to model-fitting in the Fourier transform or 
u-v plane, and fit geometrical models directly to the cal- 
ibrated visibility data by weighted least squares. 

Since the source is only marginally resolved, we must 
choose t he model a priori. We u se the same mod- 
els as in Biet enholz fc Bartell ((2005) , namely an ellip- 
tical Gaussian, which we choose to represent a pos- 
sibly elongated source such as a jet, and a spherical 
shell model a ppropriate for an optically thin super- 
nova 7 (see, e.g..lBietenholz et al.ll200 ll; IBartel et al.l2002t 
iBartel fc Bietenholz! 120081) . For the 2008 May run at 
8.4 GHz, we found that the reduced x 2 was lower if we 
excluded the data from the NL antenna, suggesting that 
there were remaining calibration problems with that an- 
tenna. Accordingly we cite the fitted sizes and uncer- 
tainties derived without using NL. (We note that the 
values obtained with the NL data are within the uncer- 
tainties, so our results do not depend on this exclusion.) 
Due to our elliptical u-v coverage, the size of the ellip- 
tical Gaussian model is more poorly constrained in an 
approximately north-south direction, and hence the ma- 
jor axis size of the elliptical Gaussian model is notably 

7 The shell model consists of the projection of a spherically sym- 
metrical, optically thin shell of emission, with the shell outer radius 
being 1.25x the inner radius. The choice of the ratio of the outer 
to inner radius has a minimal effect on our fitted outer radii. A 
ratio of outer to inner radius ~1.25 was found to be appropriate 
for SN 1993 J (Bieten holz et al.ll2003T ). and is expected on theoret- 
ical grounds {Chevalier 1982^ We further note that the 8.4 GHz 
light-curve peaked near the epoch of our first observations, sug- 
gesting that the source was just becoming optically thin. In this 
case, a disk model would be more appropriate, resulti ng in a sligh t 
underestimate of the radius (see discussion in Bartcl et al .1 120021 ). 
However the difference in angular radii is only ~8% even if the 
source is completely optically thick, and thus does not influence 
our conclusions. 



larger than the size of the circular shell model. How- 
ever, it would be a coincidence if the source were really 
elongated in the direction perpendicular to that of our 
effective resolution. 

The statistical uncertainties in the fit parameters are 
estimated from the residuals to the fit. Our final uncer- 
tainties include an additional contribution due to possi- 
ble variations in the antenna gains, estimated by artifi- 
cially varying the gains of selected antennas (e.g., Y27, 
EB, MK) and observing the resulting change in the fitted 
parameters. We discussed the uncertai nties in a similar 
fitting proc ess in more detail in e.g., iBietenholz et alJ 
(|2001f i. and IBartel et all (|2002t ). however, in the case of 
SN 2008D, the signal-to- noise ratio is relatively low, re- 
sulting in a relatively smaller contribution of the less- 
tractable systematic uncertainties. Table [2] gives the fit- 
ted angular sizes and their uncertainties. 

3.2. Determination of the Proper Motion 

We also determined the proper motion of the peak 
brightness point between our two epochs of observa- 
tions. We derive the proper motion with respect to 
the VLBA calibrator source J0919+3324, whose po- 
sition we take to be RA = 9 h 19 m 8?7871210, decl= 
33° 24' 41'.'942950 (J2000). We used the center positions 
from the fits mentioned above (which correspond closely 
to the peak brightness points in the images) as the posi- 
tion of SN 2008D. On 2008 Feb. 8, we obtained a posi- 
tion only from the 8.4 GHz observations, whereas on 2008 
May 21, we obtained positions at both 8.4 and 5.0 GHz, 
which were consistent with each other to within the 8.4- 
GHz uncertainties. 

We calculate the proper motion using the average 
of the 8.4 and 5.0 GHz positions for 2008 May 21. 
The angular displacement between our two epochs is 
+0.17 ± 0.10 and +0.17 ± 0.15 mas in RA and deck, 
respectively, where the uncertainties are intended to be 
standard errors including both the statistical uncertain- 
ties in the fi ts and a systematic contribution, estimated 
according to lPradel et a l. (2006). This displacement im- 
plies a nominal proper motion of 2.3 //as d _1 at p. a. 
45°. Since the relative uncertainty in the RA and decl. 
components of the proper motion is high, the resulting 
vector proper motion is biased and follows a Rice dis- 
tribution (assuming a Gaussian distribution for the RA 
and decl. measurements). Therefore the true proper mo- 
tion is expected to be somewhat less than the measured 
one. This is the same bias as is observed in the frac- 
tion of polarized radiation , and following the formula in 
iWardle fc Kronberd (fl97l and averaging the uncertain- 
ties in RA and decl., we estimate that the most probable 
proper motion is 2.0 //as d" 1 . We performed a Ba yesian 
analysis (see e.g.. lLore"dolll992l : IVaillancourtJl2006h treat- 
ing our anisotropic uncertainties appropriately and found 
a slightly smaller value for the most probable proper mo- 
tion of 1.6 //as d _1 , with 99.87% probability that the 
true proper motion is < 4 //as d -1 (i.e., equivalent to a 
3cr upper limit) 8 . This value corresponds to a velocity of 

8 We note that this level of positional instability could easily 
be characteristic of the reference source, J0919+2234, rather than 
due to any proper motion of SN 2008D/XRF 080109. In particu- 
lar, Bietenholz et al. (2000, 2004) found that the peak brightness 
position of M81*, the compact source in the center of M81, was 
variable on this scale over timescales of a few months. 
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TABLE 1 
VLBI Observations of SN 2008D 



Date 


Frequency 


Antennas a 


Total time b 


Bandwidth 


Recording Rate 




(GHz) 




(hr) 


MHz 


(Mbits s" 1 ) 


2008 Feb 8 


22, 8 


VLBA C 


10 


32 


256 


2008 May 21 


8, 5 


VLBA d , Ef, Y27 


12 


(51 


512 



a VLBA = ten 25 m dishes of the NRAO Very Long Baseline Array; Ef= 100 m, 
MPIfR, Effelsberg, Germany; Y27 = equivalent diameter 130 m, NRAO, near Socorro, 
NM, USA; b Maximum span in hour angle at any one antenna. HN was not available 
for this run. d PT was not available for this run. 



TABLE 2 
Angular Size of SN 2008D 



Date Frequency Spherical Shell Radius a Elliptical Gaussian FWHM b 

Angular Size 3<r Limit on /3 a pp c Angular Size 3cr Limit on /3 a pp c 
(GHz) (mas) (mas) 

2008 Feb 8 8.4 O-^tol? d 4.1 0.83jl°;|| 13 

2008 May 21 5.0 0.30^°' 35 0.79 0.98±g1S 2.5 

2008 May 21 8.4 0.0 ± 0.20 0.71 0.60 ± 0.60 2.9 



a The outer angular radius of a spherical shell model, which is appropriate for a supernova. 
The model consists of the projection of an optically thin, spherical shell with a ratio of 
outer to inner radii of 1.25. The Fourier transform of this model was fitted to the calibrated 
visibilities by least-squares. The uncertainties are standard errors (p — 84%) and include 
a contribution due to uncertainties in the antenna gains. b The FWHM major axis of an 
elliptical Gaussian model, intended to represent a jet. The ratio of the major to minor axes 
was fixed at 10, but the position angle was fitted. The Fourier transform of this model was 
fitted to the calibrated visibilities by least-squares. The uncertainties are standard errors 
and include a contribution due to uncertainties in the antenna gains. c A 3<r upper limit 
(p — 0.13% based on Gaussian statistics) on the average apparent expansion speed in units 
if c, obtained from the angular size and the age, using a distance of 28 Mpc. d At this point 
the supernova was just becoming optically thin, so would probably appear closer to a disk 
than the projection of a spherical shell. We report the radius of the spherical shell model for 
consi stency. A disk m odel would have a radius 8% larger than that of the shell model (see 
e.g., iBartel et aLll2002D . 



< 0.6 c. Conversely, if the true proper motion were zero, 
then the chance of observing one as high as we observed 
is ~15%. In short, our measurements suggest a probable 
apparent motion of ~ 0.3 c, with a value of c being very 
improbable, whereas a stationarity is compatible at the 
la level. 

4. DISCUSSION 

We have made two epochs of VLBI observations of the 
supernova 2008D, which was associated with the X-ray 
outburst XRF 080109. Our size measurements constrain 
the apparent outflow speed of any ejecta associated with 
radio emission. Our fits of a spherical shell model to 
the 2008 May data suggest a nominal average expan- 
sion speed of 54+ 57 x 10 3 km s _1 (taking the mean value 
of the measurements at 8.4 and 5.0 GHz). Supernova 
eject a usually exhibit speeds of 10 ~ 50 x 10 3 km s _1 
(e.g.. lBietenhoI3l2l)0l ), so our VLBI size measurements 
are entirely consistent with the radio emission being due 
to the interaction of a normal shell of supernova ejecta 
with the circumstellar material. Such interaction was 
found to be compatible w ith the X-ray and opt i cal ob - 
servations of SN 2008D bv lChevalier fc Franssonl (|2008fl . 
The peak ra dio brightness of SN 2008D was ~ 3 mJy 
at 8.4 GHz (|Soderberg et al.|[2008f) . corresponding to a 
spectral luminosity of ~ 3 x 10 27 erg s" 1 Hz -1 , which 
is also in the range exhibited by normal Type Ib/c su- 



pernovae (e.g., Chevalier & Fransson 2006). So the radio 
observations are fully consistent with radio emission as is 
seen from the normal interaction of the supernova ejecta 
with the CSM. 

However, as mentioned in the Introduction, it has 
been suggested t hat relativi st ic ej ecta are present 
in SN 2008D. iDado et al.l (|200l suggested that 
SN 2008D/XRF 080109 might be an ordinary, long 
GRB viewed off-axis . In this case, the models of 
iGranot fc Loebl (|2003l) predict that both the expansion 
and the apparent motion of the radio-emitting region 
should be relativistic. 

We will show that our VLBI observations provide 
strong constraints on the presence of any relativistic 
ejecta. We discuss first the expansion speed of the radio- 
emitting region. Our observations at t = 133 d firmly 
rule out an isotropic relativistic expansion, with the 3cr 
upper limit on the speed being only about ~ 0.75 c (see 
Tabled]). A GRB-like jet oriented near the line of sight 
is also ruled out, as such a jet would be expected to ap- 
pear roughly circular (e.g.. lGranotll2007fc IGranot fc Loebl 
2003) and expand with an apparent speed well over c. 

A relativistic jet with a narrow opening angle lying 
exactly in the plane of the sky would be expected to 
have an elongated morphology, and expand with an ap- 
parent speed of ~ c. A jet not exactly in the plane of 
the sky would exhibit apparently superluminal expan- 
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sion. If we conservatively assume a jet not exactly in the 
plane of the sky, of which only the approaching jet is vis- 
ible, then the major axes of the fitted elliptical Gaussian 
on 2008 May 21 suggest at best transrelativistic appar- 
ent expansion, with nominal apparent speeds of 1.2 and 
0.7 c at 5.0 and 8.4 GHz, respectively. Although our 3cr 
upper limits (/3 apP < 2.5; see Table [5]) are compatible 
with modestly relativistic expansion, even for a jet with 
a modest Lorentz factor of 5, our upper limits are only 
compatible with jets at angles > 40° to the line of sight 
because jets at smaller angles would have apparently su- 
perluminal expansion above our observational limits. In 
summary, although our measurements of the size of the 
radio-emitting region cannot conclusively rule out a rel- 
ativistic jet, they suggest more modest speeds. 

A stronger constraint on relativistic ejecta, how- 
ever, comes from our small observed proper motion. 
iGranot &; Loebl (|2003l ) show that the flux centroid of 
an off-axis GRB jet would show a large proper motion 
for periods of one month or longer. Scaling their val- 
ues to our distance of 28 Mpc, we find that the flux 
centroid would show a proper motion of > 10 /j,a,s d _1 
(corresponding to an apparent speed of > 1.6 c). Our 
3cr upper limits to the proper motion of the peak bright- 
ness point, which will be almost equivalent to the flux 
centroid, is only 4 /xas d _1 (0.6c). We can there- 
fore exclude superluminal motion between i = 30 and 
133 days. In fact our measured proper motion is con- 
sistent with SN 2008D/XRF 080109 being stationary, 
as would be expected of the emission from an ordinary, 
roughly isotropically expanding supernova, which is ex- 
pected to have a proper motion of less than a few hundred 
km s" 1 , as has been me asured, for example, for SN 1993J 
(jBietenholz et al.ll2001| ). 

We can compare our results, for example, to the 
gamma-ray burst GRB 030329, a nearby (z ~ 0.17) 
burst, whose angu lar size was measured using VLBI 
(|Tavlor et all l200i 12001 iPihlstrom et~aT1 120071 ). If it 
were in NGC 2770, GRB 030329 would have had an an- 
gular diameter of ~3.6 mas (FWHM of a Gaussian) at 
t = 83 d. This is larger than any of our 3<7 upper limits 
on the angular size. In Figure [2j we plot the size against 
the age for SN 2008D, and compare it to those measured 
for GRB 030329, and for SN 2001em, another Type Ib/c 
supernova that was a candidate for an off-axis GRB but 
was shown to have only supernova expansion velocities, 
as well as for SN 1993J, a typical Type lib supernova. 
W e note that a more re cent measurement of SN 2001cm 
bv lSchinzel et all (|2008l ) further limits the expansion ve- 
locity of SN 2001em, and an ad ditional 1.6 GHz mea- 
surement by iParagi et al.l ([2005) is consistent with the 
plotted results. Our (2a) upper limits on the extent of 
SN 2008D are considerably lower than the sizes measured 
for GRB 030329, with SN 2008D showing subluminal ex- 
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pansion, while GRB 030329 shows superluminal expan- 
sion till t > 300 d. 

The prospect of seeing an off- axis GRB event, which 
ought to be at least 10 x more common than those seen 
as GRBs, continues to prove elusive. No emission clearly 
associated with an off-axis burst has been identified, al- 
though fairly strong afterglow is emission expected from 
such j ets, particularly in t he radio (e.g.. iPerna fc Loebl 
1998j iLevinson et "alT[2002l : I Waxmanl 120041). For exam- 
ple, the Ty pe Ib/c SNe 2 002ap (|Mazzali etTai] l2002f) 
and 2003jd (jMazzali et alj 12005b were both suggested 
to possibly harbor off-axis GRB jets due to their spec- 
tral resemblance to GRB-SNe (e.g., SNe 1998bw, 2003dh 
and 20031w) In neither case, however, was the ex- 
pected late^ymfi^dto enhssion seen up to an age of ^2 
years (jSoderberg et al.|[2006f) . Indeed, a survey search- 
ing for radio emission from off-axis jets associated with 
Type Ib/c SNe found no such emission, implying that 
< 10 % of Type Ib/c SNe ar e in fact associated with GRB 
jets (jSoderberg et al.l 12006). The presence of relativistic 
ejecta was suggested both for SN 2001em and SN 2008D, 
based on their strong X-ray and/or radio emission. Both 
were bright enough in the radio to allow VLBI imaging, 
which provides the best measurements of the angular size 
and proper motion, and thus the most direct observa- 
tional constraint on the presence of relativistic ejecta. In 
both cases, the VLBI images showed no evidence for any 
radio-bright relativistic ejecta. It would seem, therefore, 
that neither high velocity optical features, nor bright ra- 
dio and X-ray emission are reliable indicators of a rela- 
tivistic GRB explosion. VLBI imaging, if it is possible, 
provides the most definitive test of the existence of a 
relativistic explosion. 

5. SUMMARY 

We report on two epochs of VLBI observations of the 
unusual supernova 2008D which was associated with the 
X-ray outburst XRF 080109. The presence of an X- 
ray outburst and the high velocities observed in early 
spectra led to the suggestion that relativistic ejecta were 
present in SN 2008D. Our VLBI observations rule out 
an isotropic relativistic outflow. Even in the case of a 
jet-like outflow near the plane of the sky, although our 
upper limit on the angular expansion is compatible with 
velocities of up to 2.5c, the upper limit on the proper 
motion of the brightness peak suggests a sub-relativistic 
velocity of < 0.6 c. 
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Fig. 1.— Image of SN 2008D at 5 GHz on 2008 May 21. The contours arc drawn at -21, 21, 30, 40, 50, 70, and 90% of the peak 
brightness of 296 £t Jy bin -1 , with the 50% contour being emphasized and the lowest contour being 3 X the image rms background brightness 
of 21 ^tjy bin -1 . The greyscale is labeled in ^tjy bm~ 1 , and the FWHM of the convolving beam (2.47 X 0.84 mas at p. a. —15°), is indicated 
at lower left. North is up and east is to the left, and the coordinate system is centered on the peak brightness, which is located at 
09 h 09 m 30?6463008, 33° 08' 20'.'124329 (J2000). Natural weighting was used to achieve the highest possible signal-to-noisc ratio in the 
image. 
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Fig. 2.— Plot of the apparent size as a function of time, comparing SN 2008D to SN 2001em, SN 1993J and GRB 030329. For SN 2008D 
and SN 2001em, we conservatively plot 2<r upper limits for the FWHM major-axis size of an elliptical Gaussian model (which are larger 
than the limits for th e spherical shell models). The u pper limits on the size of SN 2008D are from this work (see Tabled, and those for 
SN 2001em are from [Bictcnholz & Bartel (2005, 2007]). For comparison, we plot the well-deter mined expansion curve of SN 1993 J (the 
labelled dot-dashed line shows the evolution of the outer radius, taken from |Bari^l_et_al]|20 02J). The FWHM sizes for GR B 030329 are 
derived from circular Gaussian fits and an angular-size distance of 587 Mpc l lPihlstrom et al.H2007l ; ITavlor et al.| [2005. 2004), and plotted 
for their rest-frame times. For comparison, lines showing expansion with apparent speeds of 0.1 and 1.0 c area also plotted. 



